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Ground-states of spin-1 bosons in asymmetric double-wells 
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In this work we investigate the different states of a system of spin-1 bosons in two potential 
wells connected by tunneling, with spin-dependent interaction. The model utilizes the well-known 
Bose-Hubbard Hamiltonian, adding a local interaction term that depends on the modulus of the 
total spin in a well, favoring a high- or low-spin state for different signs of the coupling constant. 

We employ the concept of fidelity to detect critical values of parameters for which the ground state 
undergoes significant changes. The nature of the states is investigated through evaluation of average 
occupation numbers in the wells and of spin correlations. A more detailed analysis is done for a 
two-particle system, but a discussion of the three-particle case and some results for larger numbers 
are also presented. 

PACS numbers: 03.75.Mn, 67.85.Fg, 67.85.Hj 


I. INTRODUCTION 

The experimental realization of Bose-Einstein conden¬ 
sates (BEC) in dilute atomic gases is one of the most 
exciting recent achievements in physics, research associ¬ 
ated to this peculiar state of matter has flourished in 
late years. One remarkable development in this context 
is the realization of spinor Bose gases in optical lattices. 
In contrast to magnetic trap, where spins are frozen, in 
optical trap the atoms keep their spin degrees of free¬ 
dom. Several experimental groups have successfully cre¬ 
ated spinor BECs of 23 Na [J, @] and 87 Rb atoms (6~8|- 
Spinor gases exhibit richer quantum effects than their 
single-component counterparts, and allow to investigate 
mesoscopic magnetism. 

These experimental developments have stimulated ex¬ 
tensive study of related theoretical models mm. In par¬ 
ticular, the behavior of spin-1 bosons in a double-well po¬ 
tential can be described by a variant of the two-site Bose- 
Hubbard Hamiltonian |25| including spin-dependent in¬ 
teractions [ill! HH that affect physical properties of the 
system. Of particular interest is the case where the num¬ 
ber of atoms is small, motivated by the recent successful 
experimental trapping of few atoms with high control and 
precision [Tbi - ITsl . Remarkably, the experimental prepa¬ 
ration of only two interacting particles in a double well 
has been just reported [29j, and in principle the extension 
to threeparticles is feasible by the state-of-the art exper¬ 
iments (301 ]. These experimental achievements have gen¬ 
erated an intense theoretical effort in few-body quantum 
systems (see, for instance, [ET - fTIl ] ). Despite their sim¬ 
plicity, they still constitute a very challenging research 
field/ 

In this work we investigate the ground-state proper¬ 
ties of few spin-1 bosons in a double well. Such systems 
can be viewed as building blocks of optical lattices with 
cold atoms that can be in three different hyperfine states. 
Besides the usual Hubbard-type repulsion, the model in¬ 
cludes a spin-dependent attractive interaction between 
the particles, which may favor the establishment of a 
high or low spin state in each of the wells depending 


on the sign of the coupling constant of this interaction. 
Our basic goal is to study changes in the characteristics 
of the ground state induced by variations of the model 
parameters. 

The paper is organized as follows. In Sec. El we discuss 
the Hamiltonian and its diagonalization. A detailed anal¬ 
ysis of the two-particle system is developed in Sec. 1IIII 
where we also introduce the tools that we use to obtain 
information about relevant properties. Sec . 11VI shows rel¬ 
atively detailed results for the three-boson system, and 
discusses general trends for larger numbers, exemplifying 
with the cases of four and five particles. Final remarks 
are presented in Sec. El 


II. MODEL AND MATRIX REPRESENTATION 

Following Ref. [23|, we write a variant of the Bose- 
Hubbard Hamiltonian for a system composed of two 
wells, at positions L (left) and R (right), as 

H = e ( n L - n R ) ~ t^2(a[ a a RcT + a ] Rrj a L(J ) 

a 

+ I U Y, n i( n i -V + W 2ni )’ W 

i—L,R i=L,R 

where a\ a and a ia are the creation and annihilation op¬ 
erators of a boson in a given well (i = L,R) and in the 
spin state a = {—1, 0,1}; rii and are the number and 
total-spin operators for each well. We assume a single 
level per well in the zero-tunneling limit (t = 0), with 
energies ±e. This means that e is an asymmetry (or tilt) 
parameter, since the wells are identical only for e = 0. 
All interactions in the model are local. Apart from the 
usual (repulsive) Hubbard interaction U, the last term in 
Eq. |T|) describes a spin-dependent interaction with cou¬ 
pling constant U'. It should be noted that this term con¬ 
tributes only when more than one particles are present in 
the same well, in which case low- and high-spin states are 
favored for U' > 0 and U' < 0, respectively. The num¬ 
ber and spin operators appearing in the Hamiltonian are 
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given by 


= E 


a i<T a icr , 


Si = E« 


ha -*■ <y<y' a ia' 


( 2 ) 


where T = T x x + T y y + T 2 z, in terms of the usual spin-1 
matrices. 

From Eq. we can see that the total number of par¬ 
ticles (N t = N l + N r ), the total spin (S t = S L + S R ), as 
well as any component of the latter (which we will choose 
to be Sf) are conserved quantities. It is natural to in¬ 
vestigate the properties of a system with fixed number of 
particles. Moreover, since the energies do not depend on 
the value of S*, we can restrict our analysis to the sub¬ 
space with S* eigenvalue equal to zero, as this subspace 
is always present for any number of spin-1 particles. 

A convenient basis is provided by a set of vectors of 
the form \{N L ,N a },{S L ,S H } ,S t ), labeled by the num¬ 
ber and spin values in each well, and the total spin S t . 
Bosonic symmetry imposes that Si + Ni must be an even 
integer 23J (i = L,R). This is an interesting basis be¬ 
cause it explicitly separates subspaces of different val¬ 
ues of S t , which are not connected by the Hamiltonian. 
However, to deal with the tunneling term it is better 
to use appropriate symmetric combinations of vectors 
of type specifying 

the number of particles per spin state in each well. Since 
these two sets of vectors are related by Clebsh-Gordon 
coefficients, it is straightforward to generate any matrix 
elements. 

All results that we will discuss here were obtained 
starting from exact diagonalization of the Hamiltonian 
matrix for a given set of model parameters and parti¬ 
cle number. The on-site Coulomb repulsion U was kept 
fixed, and chosen to be the energy unit. Hence, in what 
follows we set U = 1, it being implicitly assumed that 
any quantity with dimension of energy is expressed in 
units of U. 


III. TWO-PARTICLE SYSTEM 

We begin by investigating the energy spectrum for the 
simplest case, N t = 2, with emphasis on changes in the 
ground state when varying the parameters of the Hamil¬ 
tonian. 

As an example, Fig. [T] shows the energy eigenvalues 
for two bosons as functions of the asymmetry parameter 
e, for moderately weak tunneling (t = 0.1) and spin- 
dependent coupling (U' = 0.3). For comparison, the top 
panel shows the case of decoupled wells (t = 0). We recall 
that the local energy levels are e and —e for left (L) and 
right ( R ) wells, respectively, which can be related to their 
depths. Thus, e < 0 means that the L well is deeper. 
The inversion symmetry around e = 0 is noticeable in 
the plots. 

From Fig. [1] it is clear that there are values of e near 
which the ground-state changes. These points occur at 
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FIG. 1. (Color online) Energy eigenvalues as functions of e, 
for Nt = 2 and U' = 0.3. The top panel shows decoupled wells 
(t — 0), while the bottom panel corresponds to moderately 
weak tunneling (t = 0.1). 


level crossings in the absence of tunneling, and the de¬ 
generacy lifting is stronger as tunneling increases. 


A. Ground-state changes monitored by fidelity 

Ground-state changes in finite-size systems may be 
viewed as precursors of quantum phase transitions 
(QPTs) in the macroscopic limit. Among usual tech¬ 
niques to detect QPTs, the fidelity of two ground-states 
corresponding to different sets of parameters can be used 
unambiguously for finite-size systems. This is a concept 
derived from quantum-information theory, and measures 
the similarity between two quantum states. Although 
there are generalized definitions [42j], the simplest one, 
which serves our purposes, is 

= I (*/#) I • (3) 

It defines the fidelity between any two states of the 
Hilbert space as the absolute value of their scalar prod¬ 
uct. For normalized states, 0 < T < 1. 

In our case, judging from Fig.[l| an appropriate control 
parameter is e, so that we define the fidelity 

T e (e-,N t ,t,U') = \ (e - 6 , N t ,t,U'\e + 6 , N t ,t,U') \ , (4) 

where S is small in the scale of e values, and the notation 
for the ground-state vector also includes all quantities 
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FIG. 2. (Color online) Ground-state fidelity as a function 
of the asymmetry parameter e for Nt = 2, with two sets of 
values of tunneling amplitude and spin-dependent coupling: 
t = 0.1, U' = 0.3 (solid line, right y-axis), and t = 0.005, U' = 
0.1 (dashed, left y-axis). 


that are kept fixed. Later on we will be interested in 
the effect of varying the spin-dependent interaction [/', 
in which case a more convenient choice of fidelity is 

T w (e,N t ,f,U') = | (e,N t ,t,U' — S\e,N t ,t, U' + S)\. 

(5) 

Characteristic behavior of T f for two particles is shown 
in Fig. [2] The smoother curve corresponds to the system 
parameters used in Fig. |T| Clearly defined minima are 
observed at values of the control parameter for which a 
substantial change of ground-state occurs. These minima 
become much sharper when the tunneling is reduced, as 
shown by the dashed curve, for which they appear at 
different e values because a different U' was chosen. 


B. Occupation number and spin correlations 

Although the fidelity finds parameter values for which 
the ground sate changes, it does not give direct infor¬ 
mation about the nature of states. For this we need to 
evaluate (average values of) relevant physical quantities. 
An obvious one is the occupation number of each well 
[23l |. Also interesting is the spin correlation function, 
usually employed to study magnetic properties of solids, 
as they are directly associated with magnetic susceptibil¬ 
ity. The latter not only gives the response to an applied 
field, but serves to signal the establishment of magnetic 
order. The concept of magnetic order does not make 
sense in a small system as the one that we are studying. 
However, spin correlation functions can give important 
information about the nature of the ground state with 
respect to relative orientations of the spins. 

We define the spin correlation function between the 
two wells as 

C(e, N t ,t, U') = (e, N u t, U'\S L • S fl |e, N t , t, U '), (6) 



FIG. 3. (Color online) Occupation number of the R well (solid 
line) and inter-well spin correlation function (dashed) as a 
function of e for two particles, with the same parameters that 
show sharp minima in Fig. [2] 


where and Sr denote the total spin operator associ¬ 
ated to the left and right wells, respectively. 

Figure [3] shows typical behavior of the right-well oc¬ 
cupation number and inter-well spin correlation function 
for two particles, in the same range of relative depths as 
in Fig. [2] The staircase behavior of Nr is easy to under¬ 
stand on the basis of a competition between the on-site 
repulsion, that tends to keep particles apart, and the 
energy asymmetry, that favors occupation of the deep¬ 
est well. The role of U' is to select spin states and to 
counteract the repulsion U since it is attractive for the 
appropriate spin orientations. 

By the definition ([6]), C is zero if one of the wells is 
empty or doubly occupied with zero total spin, which is 
observed in the large |e| regions of Fig. [3] In the single¬ 
occupation regime, we see that the spin correlation is 
negative, reflecting the fact that U' > 0 favors a low to¬ 
tal spin in each well, which is consistent with tunneling to 
singly occupied states with opposite spins. This is sim¬ 
ilar to the mechanism of exchange interaction between 
localized electrons of neighboring atoms in a crystal, al¬ 
though in that case, due to the exclusion principle, there 
is no need for a spin dependent interaction. 

The results shown up to now correspond to a weak 
positive U'. From Fig. [2] it is possible to see that the 
single-occupation range shrinks as U' grows. This regime 
disappears when U' exceeds a critical value U' + = 0.5, 
since the attractive effect of U' overcomes the repulsion 
U, favoring formation of a zero-spin pair in the deeper 
well. 

In the case of U' < 0, still for a two-boson system, the 
spin correlation behaves similarly to what was seen for 
U' > 0 in Fig. [3] except that the sign of C is reversed. 
Double occupancy continues to be favored by U', but 
now with maximum total spin in the doubly occupied 
well. Hence, the single-occupation region for small \U'\ 
will show ferromagnetic correlations. The critical value 
for suppression of this regime is now U'_ = —1.0 (meaning 
U'_ = -U ). 
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FIG. 4. (Color online) Regions of qualitatively different 
ground-states for two particles. The outer regions correspond 
to double occupation of the L (I,II) and R (I',II') wells. Sin¬ 
gle occupation of both wells occur in the center region, with 
AF (III) and FM (IV) spin correlations. 


C. Two-particle regimes in the (e, U') plane 


It is interesting to see a complete picture of the case 
N t = 2, varying the asymmetry parameter e and the 
spin-dependent interaction If, as shown in Fig. |4] The 
limiting lines of the various regions were obtained from 
fidelity minima (mostly F e , but the line at If = 0 is bet¬ 
ter seen with Fit). I 11 regions I and II there is double 
occupation of the L well, with Sl = 0 and 2, respec¬ 
tively. Regions I' and II' are equivalent to I and II, but 
the particles are located in the R side. States of one par¬ 
ticle in each well are observed in the two central regions, 
III and IV, respectively with negative and positive spin 
correlations. Borrowing denominations from magnetism, 
we can say that spin correlations are antiferromagnetic 
(AF) in region III and ferromagnetic (FM) in IV. 

It is worth observing that the characterization of states 
in terms of doubly or singly occupied wells refers to the 
average occupation. Fluctuations occur for any t ^ 0, so 
that the sharp boundary lines of Fig. |4] are only sharp 
for very weak tunneling. For instance, choosing a fa¬ 
vorable condition for single occupancy inside region III 
(e.g, symmetric wells and small \U'\), a possible mea¬ 
sure of the probability of double occupancy is given by 
T> = 1 — Jy 2 , where Ft is the fidelity between the ground- 
states with zero and nonzero tunneling. This is shown 
as a function of U/t in Fig. [5l where we can see that it 
is null for decoupled wells, and grows with the tunneling 
strength, approaching 1/2 for large t, when all possible 
occupancies become energetically equivalent. The bot¬ 
tom curve ( U ' = 0) is very similar to the one obtained 
for two fermions in a double-well (2fj. Here, the effect of 
increasing U' is to lower the energy of doubly occupied 
wells. 



FIG. 5. (Color online) Variation of the average double¬ 
occupancy as a function of the the ratio between Hubbard 
repulsion and tunneling strength in the two-particle ground- 
state for a symmetric double-well. The curves correspond to 
values of the spin-dependent interaction within the region III 
of Fig. |31 



FIG. 6. (Color online) Fidelity (continuous) and A-well oc¬ 
cupation (dashed) for Nt = 3, t = 0.005, and U' = 0.1 


IV. THREE OR MORE PARTICLES 

Adding a third particle to the system leads to notice¬ 
able qualitative changes. Figure |6] shows the fidelity and 
R-well occupation for varying e. The plot is for U' > 0, 
but the behavior for negative U' is qualitatively the same. 
Notice that there is a central minimum of the fidelity, 
corresponding to the extra step observed for Nr, but no 
central plateau, in contrast with Fig. [2] since it is im¬ 
possible to have the same number of particles in both 
wells. The most striking difference comes from the spin 
correlation function, shown in Fig. |T| For positive U' 
the correlation is negative but close to zero except, for a 
slight increase in magnitude near the changes of ground- 
state. This is due to the fact that there is always either 
zero or two particles in one of the wells, which means 
that one of the spins is always essentially null. In con¬ 
trast, for If < 0 the whole central region, where none of 
the wells is empty, presents a large positive correlation 
between the spins of the two wells. 

As in the two-particle case, when the magnitude of a 
negative If grows beyond a critical value (once more, 
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FIG. 7. (Color online) Spin correlation functions for Nt = 3 
and t = 0.005 in the cases of positive and negative U'. Notice 
the enlarged scale and displaced origin for the solid line. 



FIG. 8. (Color online) Variation of the ground-state fidelity 
with t and U’ for t = 0.1 in a three-particle system. Notice 
that the fidelity axis is inverted, and the values are normalized 
to fall in the range [0,1]. 


U'_ = —1), all three particles stick together in one of the 
wells (with maximum spin S) = 3), so that the two in¬ 
ner regions around the central minimum in Fig. [G] cease 
to exist. On the other hand, for positive U' the fidelity 
minima are located at the same e values seen in Fig. [Gj 
independently of U 1 . This is shown in Fig. [51 where we 
plot a normalized fidelity T n (rescaled to fall in the range 
[0,1]) as a function of both e and U'. The observed inde¬ 
pendence on V is due to the constraint that N + Si must 
be even, which means that the minimum total spin in a 
given well is Si = 0 for Ni =2 and Sj = 1 for N = 3. It 
then follows that the eigenvalue associated with the U' 
term [see Eq. ©] is Si(Si + 1) — 2 Ni = —4 in both cases, 
so that the change in ground-state is driven by a balance 
between e and the local repulsion U. 

A ground-state diagram like that of Fig.[4]can be built 
for N = 3, as shown in Fig. [9] Here too the regions 
labeled with primed roman numbers are equivalent to the 
corresponding unprimed ones upon the exchange L tA R. 



FIG. 9. (Color online) Regions of qualitatively different 
ground-states for three particles. The outer regions corre¬ 
spond to double occupation of the L (I,II) and R (I',II') wells. 
In the central regions (III, IlF, IV, IV*) there is single occu¬ 
pation in one well and double in the other, with very weak 
AF correlations for U' > 0 and strong FM correlations for 
U' < 0. 


In I and II we have essentially three particles in L while R 
is nearly empty, with St = 3 in I and St = 1 in II. Regions 
III and IV correspond to double occupation of L and the 
third particle in R, with the spin states compatible with 
spin correlations as shown in Fig. |T| 


A. Larger N 

The main trends in behavior observed for two and three 
particles appear in general for even and odd particle num¬ 
ber. This is exemplified in Fig. [TU] for N t = 4 and 5. Ob¬ 
viously, the staircase increase in single-well average occu¬ 
pation has a number of steps reflecting the total particle 
number. Each step is accompanied by a sharp minimum 
of the fidelity. 

Spin correlation functions are shown in Fig. [TO] only 
for U' > 0. They reproduce the scenarios already viewed 
for two and three particles. Also similarly to those cases, 
correlations for U' < 0 are positive (when not null), and 
have significant values, that remain nearly constant be¬ 
tween jumps of single-well occupations. 

The “phase diagrams” are similar to Fig. [4] for even 
N t and to Fig.[9]for odd N t , except that the central part 
has a growing number of nested regions as the number of 
possible distributed occupations increases. 


V. CONCLUSIONS 

We studied the possible ground states of spin-1 bosons 
in double-well potentials. Such systems can model the 
basic unit of optical lattices with trapped cold atoms, for 
which parameters such as the depth of the wells, ampli¬ 
tude of tunneling between them and interactions between 
particles can be controlled. We based our analysis on the 
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FIG. 10. (Color online) Normalized fidelity (J- n ), spin correlation function (C), and A-well occupation for systems of four 
particles (left panel) and five particles (right panel), with t = 0.005 e U' = 0.1. 


usual Bose-Hubbard Hamiltonian with an additional on¬ 
site spin-dependent interaction, as proposed in Ref. [23( |. 

Even considering different total numbers of particles, 
restriction to two wells allowed us to exactly diagonal¬ 
ize the Hamiltonian matrix in relevant subspaces, and to 
study changes of the ground state induced by variation of 
the model parameters. Our analysis focused primarily on 
a regime of weak tunneling relative to the local repulsive 
(Hubbard) interaction, which we kept fixed. The variable 
parameters were then the relative depth between the two 
wells and the local spin-dependent coupling, which we 
allowed to be positive or negative, respectively favoring 
low and high total spin in each well. 

We showed that regime changes in the system can be 
detected by evaluating the ground-state fidelity as some 
parameter is varied. This quantity presents sharp min¬ 
ima at parameter values for which the nature of the 
ground state changes, which should correspond to critical 
values for occurrence of quantum phase transitions in the 
macroscopic limit. Information on the nature of the dif¬ 
ferent ground states has to be sought through evaluation 
of average values of appropriate physical quantities, like 
the number of particles and total spin in a given well, or 
inter-well spin correlations. 

The results allowed us to construct maps of different 
regimes in parameter space, identifying regions of full 
occupancy of a single well, and regions with particles 
distributed in both wells, in which case either ferromag¬ 


netic or antiferromagnetic spin correlations between the 
two wells occur as a consequence of the spin-dependent 
interaction. This was done in detail for system of two 
and three particles. Comparison with some results for 
four and five particles revealed that the main qualitative 
differences occur between even and odd total number of 
particles in the system. 

This study can be extended to larger numbers, not 
only of particles but also of wells, the limitations being 
only computational. Work in dynamical processes in the 
same model is now in progress, including investigation of 
spin effects on transistor-like behavior |43[ in the three- 
well case. Additionally, by increasing the number of wells 
we can address the problem of Anderson localization 01 
with a distribution of well depths, taking into account the 
effect of nonzero spin and spin-dependent interactions. 
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